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ABSTRACT 

Context. Cool star model atmospheres are a common tool for the investigation of stellar masses, ages and elemental abundance 
composition. Theoretical atmospheric models strongly depend on the atomic data used when calculating them. 
Aims. We present the changes in flux and temperature stratification when changing from iron data computed by R.L. Kurucz in the 
mid 90s to the Kurucz 2009 iron computations. 

Methods. MAFAGS-OS opacity sampling atmospheres were recomputed with Kurucz 2009 iron atomic data as implemented in the 
VALD database by Ryabchikova. Temperature stratification and emergent flux distribution of the new version, called MAFAGS-OS9, 
is compared to the former version and to solar flux measurements. 

Results. Using the Kurucz line lists converted into the VALD format and new bound-free opacities for Mg i and Al i leads to changes 
in the solar temperature stratification by not more than 28 K. At the same time, the calculated solar flux distribution shows significantly 
better agreement between observations and theoretical solar models. 

These changes in the temperature stratification of the corresponding models are small, but nevertheless of a magnitude that affects 
stellar parameter determinations and abundance analysis. 

Key words. Atomic data - Sun: abundances - Stars: abundances - Stars: atmospheres - Stars: fundamental parameters - Stars: 
late-type 



^ ■ 1. Introduction IMashonkina et al J d2008bl). IMashonkina et al J d2008al) . IShi etaH 

d2008l) and lMashonkina et all d2009l) . 

Theoretical model atmospheres are the most basic tool in the in- For basic work on stellar parameters, that also enters the 

Cj ■ vestigation of stellar parameters such as temperature, mass, age abundance analy sis via the selection of th e proper mode l, 

and of the element abundances present in stellar atmospheres. Megessie i] d 1 998b ICastelli & Kuruczl (11994. ICastellil d 19991) , 

Thus the accuracy in determining all these stellar properties Grupp (2004bT) andlGrupp & Mashonkinal d2008h show the influ- 

strongly depends on the accuracy and reliability of the under- ence of different model atmospheres on various stellar parameter 

lying stellar atmosphere. determination methods. 

Many studies based on model atmospheres have shown the As a basic ingredient, bound-bound and bound-free atomic 
usefulness of model atmospheres for stellar paramete r and abun- processes dominate the absorption properties of late type stellar 
dance work. The work of Fuhrmann and col l eagues ( Fuhrmann atmospheres over a wide spectral range. Taking into account as 
d 1 9981). IFuhrmannl d 1 9991) . IFuhrmannl (2002), Fuhrmann (2004) complete as possible datasets for these major opacity sources is 
and IFuhrmannl d2008l) ) on nearby stars of the galactic disk shown to be vital for accurate stellar parameter and abundance 
and halo might serve as an example of the validity of this analysis work, 
approach. While Fuhrmann stays with the assumption of lo- 
cal thermal equilibrium (LTE), many precision studies of el- 
ement abundances in cool stars remove this simplification 2. The model atmosphere code 
and calculate element abundances based on full radiative 

and collisional equilibria (Non-LTE). For examples, see e.g. Our study is based on the MAFAGS-OS theoretical model atmo- 

sphere code described by Grupp (2004a). Based o n the PDF ver - 
sion of T. Gehren in the reprogrammed version o f lReildd 19871) , 
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Plane-parallel ID geometry. 

Chemical homogeneity throughout the atmosphere. 
Hydrostatic equilibrium. 

Convection is treat ed acco rding to the formalism of 
ICanuto & Mazzitellil d 19911) and ICanuto & Mazzitellil 
19921) . 



- No chromosphere or corona. 

- Local thermal equilibrium. 

- Flux conservation throughout all 80 layers. 

While these assumptions might break down for hot stars, stars 
with very extended atmospheres and the coolest stars, they are 
valid for stars in the range of T e ff « 4000 • • • 15000K and for 
gravities from the main sequence down to log(g) as 0. 

3. Atomic data 

The main target of this study is to investigate the effect of the new 
Kurucz iron line data on models and energy distributions. These 
data were com piled for us in a p rivate data extraction run us- 
ing the VALD ( Kupkaet al.ll 19991) database in January 2009 by 
Tanya Ryabchikova. This private extraction run was performed 
including predicted levels in all elements of the iron group. 

3. 1. Old and new line list 



8x10 



1.5X10 7 - 1 " 6x10 



4x10 



(!) 

E 



E 



.0x10 



7 z 2x10 




a) 



bnistation state 



5.0x1 fa 



1 1 1 r r i 



20 22 24 26 28 

Element number 



the levels greatly increased. The line list are completed for the 
rest of t he iron group u s ing ob serve d and predicted levels of 
iKuruczj d 1 994al) . Ikuruczl d 1 994d) and lKuruczl dl994bl) . Both the 
new data of iron and the data for the rest of the iron group el- 
ements from the stated Kurucz CD-roms is incorporated in the 
VALD database used for data extraction by Tanya Ryabchikova. 
All lin es were processed using the data extraction algorithms of 
VALD dKupka et al.lll999h . 

The s e data replace ol der d ata comp i led in 1998 from lKuruczl 
dl994al) . IKuruczj d!994d) . and lKuruczl d!994bl) . For the old line 
lists lines between both observed and predicted levels have been 
extract making use of the files GFxxyy.GAM and GFxxyy.LIN 
where xx stands for the atomic numbers 20 through 28 and yy 
codes the ionisation states 00, 01 and 02. We used our own re- 
duction software adopted from the software delivered with the 
CDs to process these files. 



For the non iron group elements, Kurucz & Bell ( 1995) data 
were used in the old line list. iKurucz & Belli (|l995) data was 
also included in the new list, but also its continuous updates in- 
corporated in VAL D. For detailed references to the new list see 
iKupka et al.l (11999) and the compilations on the VALD Internet 
platforrry. Both lists - old and new - contain diatomic molecules 
taken from lKuruczl dl993l) and lKuruczl d 19991) . 

FigureQ]a) shows the number of lines in the final list for the 
iron group elements. The enormous increase of iron lines is ob- 
vious. Small deviations in the number of Cr lines are due to the 
fact that in the old list, Cr lines had been included both below 
500Aand above 100000A. These regions are unimportant for the 
solar line blending and the bug has been corrected in the new 
version of the line list. 
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Fig. 2. Number of lines for the elements 2 ■ • ■ 52 on logarithmic 
scale. Color coding is the same as in fig.Q] 



Fig. 1. Number of lines for iron group elements in new (black) 
and old (grey) MAFAGS-OS line lists. 

a) For the elements Ca through Ni, sum of ionization states I,U 
and III. 

b) For the first three ionization states of Fe. 



Our new line list was compiled based on the abovemen- 
tioned dat a extraction of V ALD and includes the latest Fe i to 
in data of Kurucz (2009), taking into account observed and 
predicted transitions for each ion. These new calculations for 
iron included additional laboratory levels and more configura- 
tions than earlier work. Thus, the number of lines connecting 



The number of lines for the non iron group elements Li 
through Te is shown in fig. [2] It becomes obvious that there are 
few changes for all elements but iron between both lists. 

3.2. Changes in continuous opacities 

Another change of atomic dat a with respect to the 20 04 model 
is the inclusion of TOPBASE dCunto & Mendozall 19921) bound- 
free cross section s for Al i dMendoza et alj 1 19951) and Mgi 
(But ler et alJ [T990). Old and new models make use of Bautista 
dl997n " ionization cross sections for Fei (See fig.|4|. 
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3.3. Line selection 

Our model considers the ionization states one, two and three 
and runs with m 86000 wavelength points between 500A and 
850000A. TableQ] shows the number of bound-bound transitions 
in the new and old line-list together with the number of lines that 
were selected by the code to sample our standard solar model 
opaciti es (For a description of the selection criteria see [Grupp 
2004a). The number of selected lines corresponds to our stan- 
dard solar model. Both line lists contain additional 4247857 lines 
related to transitions in diatomic molecules. 

Table 1. Data for the new and old atomic linelist. 



List 


Total number 


Number of selected 




of atomic lines 


lines for sampling 


1998 


16012015 


322040 


2009 


35352145 


621406 



4. Solar flux and temperature structure 

Figure|3] shows the comparison o f the old and new m odels 
with observed solar flux data from iNeckel & Labsl ( 1 19841) and 
IWoods et al.l (Q996). Despite remaining inconsistencies in the 
wavelength range blueward of 5000A, it is obvious that the 
model computed using the new iron data shows both less scatter 
and better absolute agreement with observed quantities. 

Nevertheless, the region between 3000A and 5000A shows 
remaining shortcomings of about +10% in flux, that reach +20% 
at some points. These discrepancies and the slight underestima- 
tion of the solar flux in the very scatter-dominated region below 
3000A are present in both datasets. Still, the computations with 
the new iron data are closer to observation even in those regions. 

The flux changes are mai nly due to the new da ta for iron. 
The influence of new All (Mendoz a et al.l fl995) and Mgi 
dButler et al.l fl990) bound-free cross sections is shown in Fig. [4] 
The flux differences are small in the visible wavelength range 
and contribute only to some larger deviations in the near ultra- 
violet region. These changes in the UV do not significantly af- 
fect the atmosphere stratification, but they can be important for 
non-LTE calculations, as the departure coefficients of levels with 
transitions in the UV might well depend stro ngly on photon pro- 
cesses in this region (see e.g. lShi et al.l d2008)). The rest of the in- 
put atomic data concerning bound-free transitions is unchanged 
with respect to the "old" model. 

Figure[5] shows th e comparison of our la test model to 
the ATLAS9 model of ICastelli & Kuruczl d2003l) . The selected 
ATLAS model is the ASUNODFNEW modefl computed with 
new ODFs and the overshooting option turned off. 

The ATLAS model chosen uses the new ODF tables of 
Kurucz and has the overshooting option turned off. While both 
models show similar flux structure in the red and near IR region, 
with the ATLAS model having slightly less flux here, they sig- 
nificantly differ in the UV to green part. This leads to the ATLAS 
model underestimating the observed solar flux in the red and IR 
by as 5% at 9000 A and to a good agreement between ATLAS 
and measured solar flux in the green spectral range. Looking at 
the picture as a whole, neither of the two versions fits the solar 
flux significantly better than the other, but both have good and 
not so good spectral regions. The difference between the models 

2 http://wwwuser.oat.ts.astro.it/castelli/grids.html 
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Fig. 4. Comparison of MAFAGS-OS model atmospheres using 
"old" boun d free data for Al i and Mg i and recent detaile d cal- 
culations of lMendoza et al.l (Il995l) and lButler etaD (Il990l) . 

is quite likely to emerge from the lack of the new iron lines in the 
ATLAS9-NEWODFs and some different UV continuous opacity 
sources used in both models. In addition, the models use differ- 
ent formalisms to treat convection, which will also influence the 
overall flux distribution. 

The changes in solar temperature structure due to the two 
different line lists are shown in Fig. [6] The overall effect is small, 
nevertheless changes of « 5 ■ • ■ 25 K in the line forming regions 
will slightly affect line formation. 

5. Discussion and conclusions 

We have shown the influence of the new iron Kurucz data on 
the solar model and solar energy distribution. This new data 
shows significantly better agreement with the observed solar flux 
with respect to scatter and absolute values. Changing to the new 
data also slightly changes the solar temperature stratification. 
Therefore, the new data are shown to be a real improvement of 
the model and helps to better represent solar atmospheric ob- 
servables. 

Remaining discrepancies are mainly found to be in the short 
wavelength part below the solar flux maximum. They might 
be related to missing or inaccurate data in the resonances of 
strong bound-free absorbers or slight misplacements of the cen- 
ter wavelengths of strong lines. 

Some major implications of the improved representation of 
the stellar flux are: 

- Changes in synthetic colors and line indices. 

- Changes in abundance analysis due to changes in tempera- 
ture stratification of the model atmosphere. 
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Fig. 3. Comparison of observed (black) and model solar flux for the old (blue) and new (red) solar models. Observations are taken 
from iNeckel & Labsl (ll984l) for the spectral range above 3300Aand from Neck el & Labsl d 19841) for the near UV region below 
3300A. 

Top panel: Direct comparison of solar flux data. 

Mid panel: Deviation of the models from the measured solar flux in absolute values. 
Bottom panel: Deviation of the models from the measured solar flux in relative values. 



- Changes in non-LTE analysis due to changes in the UV 
fluxes that determine photo-ionization processes emerging 
from transitions in this spectral region. 

- Changes in stellar parameters if they rely on synthetic colors, 
line indices or individual lines. 

- The lower scatter also allows better automatic comparison 
of synthetic and observed spectral data. A work field that 
shows incr easing importanc e with large survey project s such 
as SPS S (lYork et al'll2000h LAMOST dSu et al.ll 19981) and 
GAIA dPerrvman et al.ll200ll) on the way. For those projects 
the mere number of spectra acquired does not allow for man- 
ual inspection, requiring automated procedures of spectral 
analysis. 

Although the impact on the temperature stratification is rela- 
tively small, changing to the new data significantly changes the 



flux of the solar model, placing it into better agreement with ob- 
servations. 
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Fig. 5. Comparison of the new MAFAGS-OS model (red) solar 

atmosphere with the recent ATLAS model with NewODF and 

the overshooting option switched off (blue). 

Top panel: Direct comparison of solar flux data. 

Mid panel: Deviation of the ATLAS model from the MAFAGS- 

OS9 model in absolute values. 

Bottom panel: Deviation of the ATLAS model from the 
MAFAGS-OS9 model in relative values. 
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